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Abstract 

This work is a theoretical investigation on the spin-polariton (polarized single photon) entangle- 
ment in nitrogen vacancy (NV) centers in diamond in order to interpret the results of two landmark 
experiments |4l |60] published in Science and Nature. A Jaynes-Cummings model is applied to analyze 
the off- and on-resonant dynamics of the electronic spin and polarized photon system. Combined 
with the analysis on the NV center's electron structure and transition rules, this model consistently 
explained the Faraday effect, Optical Stark effect, pulse echo technology and energy level engineering 
technology in the way to realize the spin-polariton entanglement in diamond. All theoretical results 
are consistent well with the reported phenomena and data. 

This essay essentially aims at applying the fundamental skills the author has learned in Quantum 
Optics and Nonlinear Optics, especially to the interesting materials not covered in class, in assignments 
and examinations, such as calculations on matrix form of Hamiltonian, quantum optical dynamics with 
dressed state analysis, entanglement and so on. 

1 Introduction 

1.1 Investigation on entanglement towards photonic applications 

At the very heart of applications such as quantum cryptography, computation and teleportation lies 
a fascinating phenomenon known as "entanglement" -the spooky, distance-defying link that can form 
between objects such as atoms even when they are completely shielded from one another [8]. This type of 
correlation between particles is the deepest difference between Quantum and Classical world [TS] [TU] [7], 
and is also the key to realize "Qubits" , the unit of Quantum Information and future Quantum Computer. 
Because the photon is the best particle to carry information and to propagate it to a distant receiver at the 
fastest speed, realizing photon-photon and multi-photon entanglement is the best choice for future photonic 
applications [29 . Now scientists have realized two-, four-, six-, and even more photon entanglement using 
parametric down-conversion and other nonlinear optical technologies (see for example 149| ). However, 
it is hard to operate on photons directly, so it is necessary to study other particles' entanglement and 
their coupling with photons. Until now, successful quantum entanglement has also been demonstrated 
with individual ions and atoms (see for example this review paper [SD]), electrons in superconductors 
[1], and coupling superconductor qubits to cavity light mode [38 . However all of these apparatuses are 
either too big for application or dependent on tough conditions (like ultralow temperature) to maintain 
the entangled states, in other words , they are far from practical applications. Only in recent years have 
people made considerable progress in acceptable conditions with integrable entangled solid state systems 
[2SlllO]~like semiconductor Quantum dots [44i, photonic crystals |67 -and has "spintronics" come into the 
science and technology community [24] [55] . 

It turns out electronic spin is a good carrier that can be used for computing and storing information 
at the same time [S5] [SS] [3] [7D]. What's more, it is also easy to be coupled to polarized photons, which 
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is believed to be a better way to enlarge our information transmitting capability without costing more 
energy and fibers [6T] [35] [45]. For the practical application in quantum memory for future computer, 
it both asks for long dephasing or operating time (T2) and acceptable working frequency, as well as the 
capability of operating in room temperature. These properties are the materials' intrinsic properties, and 
will hardly be modulated along with the development of science and technology. A comparison among 
several promising candidates for quantum computer and quantum information applications is listed in 
TableO 

Table 1 ' Current performance of various qubits 
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One qubit 
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Infrared photon 
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0.47* 
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spins bound to '"P:^*Si 


0.6 s 


S 




49 


^Si nuclear spins in ^*Si 


25 s 


5 




50 


NV centre in diamond 


2 ms 


2 
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60, 61, 65 


Superconducttrtg circuit 


4|is 


o./t 


10* 


73, 79, 81, 109 



Measured 77tinnesarestown,eJcc^ for photons where 7? is replaced by twice ttwhold-ttrnc 
(compa fable to ?i) of a letecomrTnunicalkjn. wavelength photon ifi fibre. Benchmarkirig values 
shgw approximate error rates far single or multi-qubitgates. Value?! marked wit ti asterisks are 
found by quanlufii process or slate tomography, ar>d give the departure of the fidelity frorn 
100%. Values marked with daggers are found with randomijedbenchmarking"*'.Ot tier values 
are rough expefimental gate error estimates. In the case ot photon^ Iwo-qubit gates fail 
frequently but success is heraldect error raEes shown are conditional on a fierakled Success. 
NV, nitrogen vacancy. 

Figure 1: Current performance of various qubits. Reprinted by permission from Macmillan Publishers 
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As shown in the table NV (nitrogen vacancy) center in diamond has almost the longest operating 
time slot and can work from microwave to visual light band. It also can run at room temperature, 
owing to the Zero Phonon Line emitting ^54j . The probability of building a quantum computer basing on 
diamond and defects, like the NV center, is well confirmed by many research groups [251 15111 1551 [5S] . 
In recent years, the electron structure of NV center in diamond is widely studied using different methods 
[111 [ig 1211 [221 [23 Ell |34l ISa [eg , and the difficulties of entanglements among the spaced NV centers [5] 
and coupling to superconductor circuits [5] , integrated to photonic cavities [SHI US and microrings [IH] 
are also overcome very recently. 

However, it is a real challenge to realize the entanglement between electronic spin and polarized photon 
in NV centers in diamond [5 2) [2]. There are several difficulties. For example, the energy bands are very 
complex compared with individual ions and atoms exposing to the crystal environment. And it is hard to 
modulate the energy level structures to include two electronic spin states with almost the same energy gap 
to transfer to the common ground or excited states to easily operate and generate entanglement between 
photons and spins. In 2010, coherent operation and entanglement between electronic spins in NV center 
and polarized photons are performed successfully by Buckley and Togan and their teams ^ |60j . 

Emre Togan and collaborators reported their work in August 2010 [6^, first-time showed the stable 
entanglement between spin state and photon polarization in a solid-state material which potentially can 
be demonstrated in room temperature. Buckley's work was reported in Nov 2010 [4J, on coherence control 
experiment for single spin in diamond, with the Faraday Effect (FE) and the Optical Stark Effect (OSE) 
observed under scheduled exciting pulse series. FE is usually a tool to control polarization rotating of 
emitted light from solid system. And OSE can cause energy band structure adjust and luminescence 
spectrum shift if an ultrafast (usually in femtosecond scale) laser pulse is applied on the sample. The 
OSE has been investigated both theoretically and experimentally in QDs of semiconductor systems (see. 
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for example, the first experimental observation of OSE in semiconductor QD 62 ) and other solid optical 
systems [36l [57l |65] , and has been used as a tool to control single electron and photon or phonon interaction 
in solid nanostructures in recent years |321 133j . These two effects are important for quantum operations 
(see for example [S^ [HI)- 

In Togan and Buckley's work, they all used the A2 singlet state as the excited state and a pair of 
lower rris = ±1 magnetic momentum states as the ground states, and operated the system with tunable 
637nm lasers. Their results show a promising way toward diamond, or more generally, solid-system based, 
quantum logical gates, multi-spin-photon entangled network, and other practical applications in related 
fields. This article will mainly focus on these two experimental reports and represent their results using 
a unified model. 

1.2 A brief introduction of entangled states 

Before we move onto the in-depth discussion, let's briefly introduce the basic theory of entanglement in 
general. 

As discussed in Gerry and Knight's book [20], the four Bell states are given by 



|$+; 
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where "H" and "V" are orthogonal states with "opposite" properties, for example, spin-up versus spin- 
down, polarized-to-z versus polarized-perpendicular-to-z, and so on. If a quantum system is in one of the 
four Bell states, then we can say this system is in a "maximally entangled" state. "Maximally entangled" 
means that when we trace over quantum substate B to find the density operator pA of quantum substate 
A, we obtain a multiple of the identity operator. For example, let's consider a polarized photon A and a 
spinning electron B are in the state of l^*^), we have 

PA^trB{\^+){^+\)^ J2 I*+>(*+N)bb(*| 

= ^mAA{H\ + \V)AA{V\) = llA, 

similarly, pB = This means that if we measure photon A polarized along any axis, the result 

is completely random, we find polarization parallel to the axis with probability 1/2 and polarization 
perpendicular to the axis with probability 1/2. This is also true for measuring B's spin state. Therefore, 
if we perform any local measurement of A or B, we acquire no information about the preparation of the 
state, instead we merely generate a random bit (number set combined with or 1). 

However, when we repeat our measurements on A and B, if we get the state of A, then the state of B 
is acquired. Because A and B have some correlation. As Gerry and Knight's definition, the correlation 



function can be written as 

C{H, V) = Average[A{H)B{V)], (3) 

which means the average of probing A in state H while B in state V. Here, A is in H then A{H) — 1, 
otherwise, A is not in H, then A{H) — —1. Similar to B(V). For the 1^'"'') state's case, we can write 

C{H,V) = Pr{\H)A\V)B) + Pri\V)A\H)B) - Pr{\H)A\H)B) - Pr{\V) a\V) b), (4) 
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and 

Pr{\H)A\V)B) =a{H\b{V\^+) = 1/2, 
Pr{\V)A\H)B) = 1/2, 
Pr(|lf)^|F)B) = 0, 
Pri\V)A\V)B)^0, 

thus that C(i7, 1^) = 1. That means when we know photon A is in H state, then we are sure that spin 
B is in state V. Because H and V is different, we call this entanglement as non-parity or antialigned 
entanglement. In contrast, there are parity or aligned entanglement, if A and B are always with the same 
type of states. This conditional probability can be measured and verified among the statistics on a large 
number sets of measurement experiments. If we span our measurement basis in the form of Bell states, we 
can get the other correlation function for conditional measurement. And if any of the correlation function 
gives a value greater than 0.5, then the system is in an entangled state. 

For a full description of entanglement, entropy and negativity are important conceptions, you can find 
more details in [64]. [12] and [48] give a systematical introduction on multi-particle (bipartite, cluster, 
GHZ and so on) entanglement and measurement methods. And besides the discrete type (like spin states) 
entanglement, some continuous variables can also be entangled, like momentum, position, energy and so 
on. A theoretical introduction can be found in [T3] [S3] [T^ . 

Our discussion in this essay is basically on the theory and experimental methods to establish a |^+) 
type entanglement between a polariton and an electronic spin in negatively charged NV center of diamond. 
The essay, therefore, includes the following topics: a simple but unified Quantum Optical theory to 
describe the polariton-spin interaction (both off-resonant and on-resonant dynamics) in NV center of 
diamond, the NV center's electronic structure and possible level structure to perform a good spin-polariton 
entanglement, represent one recent experiment-its methods and results-using our model. Through the 
course of theoretical analysis, I will also explain the phenomena of Faraday effect (FE) and Optical Stark 
effect (OSE), which are inevitable effects to consider, and point out some possible errors in two highly 
impacted articles. 

2 Quantum Optical theory of spin-light interaction 

From this part we mainly reference and compared with the experimental results and theoretical analysis 
of two up-to-date articles and one perspective article: [4], pubhshed in Nov 2010, [60], published in 
Aug 2010, and [35] as a comment on [3], published in Science in Nov 2010. I will build up a unified 
theory frame mainly inspired by Buckley's article, using the knowledge I have learned through the course, 
to explain the reported phenomena of Faraday effect (FE), Optical Stark effect (OSE), spin echo and 
entanglement generation between nitrogen vacancy center electronic spin and photon, and to point out 
some new discoveries and argument through comparing my calculation with published results. To make 
our discussion well readable, some contents are cited from the original papers without notice, and some 
parts will be highlighted with colored fonts to distinguish my disparate arguments with the publication. 

2.1 Off-resonant dynamics theory and FE, OSE in diamond 

Now let's consider the nitrogen- vacancy center under a coherent laser pump and excited from ground state 
to excited state. The initial and final states for the joint spin-photon system can be written as 

\A"''^)^\9,)\n + l) (6a) 
|4"'^)) = |e,)|n), (6b) 

where \gj)i\ej)) are the bare ground (excited) states of the NV center orbital transition to spin number 
iTT-s — j, j — ~li 0, -1-1 state, and \n) is a photon-number state of the electromagnetic field, n = 0, 1, 2, • • • . 
As the energy splitting between different rus states is so small compared with the detuning, we can safely 
simplify the interested energy levels into a two-level system, with a formalism similar to Jaynes-Cummings 
model. The energy level structure is shown in Fig[2] 
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Figure 2: Level structure for spin-photon interaction in diamond using Jaynes-Cummings model [4j. A 
shows the simplified levels and energy shift with light-electron interaction. B,C show the detailed level 
structure in reality. From B. Buckley, G. Fuchs, L. Bassett and D. Awschalom. Spin-Light Cohe rence] 
for Single-Spin Me asurement and Control in Diamond Science, 330 (2010), 1212[ 1]. Reprinted with 
permission from A A AS. 



Now we can try using dipole approximation to write the interaction Hamiltonian as Hint 

-(a^i7_ — ad+), (7b) 



= CflE (7a) 



where I used fl — \ ((7+ + CT-) and E = i\fl\ {a^ + a) for the electric field and optical dipole, respectively, 
and a^(a) and (T+((7_) are creation (annihilation) operators for optical photons and NV center excitations, 
respectively, and we neglect energy-nonconserving terms a(T_,a^d'4. as a result of rotating wave approx- 
imation. For this case, I also make C = F^-^^i Fdw = 0.04 ± 0.01 is the Debye- Waller factor which 
empirically accounts for the reduced resonant coupling between NV center ground and excited states due 
to non-resonant phonon-assisted transitions. And 

where Ep^ = 1.945 is the photon energy and tid = 2.4 is the refractive index of diamond. We also 
used the on-resonance optical Rabi frequency as 

^^0 = ^^^ImPo|cos(,S), (9) 

where cos(5) accounts for the geometric coupling between the NV center dipole and the linearly polarized 
light. 

Considering the non-interacting Hamiltonian for the spin and light field, given by 

Ho = Eph{a^a+]^) + E,^-, (10) 

where Ej is the transition energy for the spin state with spin number nis = j, = \ej){ej \ — \gj){gj\ and 
we have ignored the zero-field energy and make the average spin energy to be zero for simplicity, we can 
finally write the Jaynes-Cummings-like Hamiltonian for the system as 

Hjc^Ha + H,nt (11a) 
= Epha^lx + Ej'^ + ^{da+ + a)&_). (lib) 
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Now let us rewrite the Hamiltonian into matrix formalism. Using the following relationships: 

0'+ = |ei)(ffj|, (12a) 

a-^\g,){e,\, (12b) 

^z^\ej){ej\-\g.j){gj\, (12c) 

a\n) = ^/Ti\n - I) , (12d) 

a)\n) = V?i + l|?i + 1) , (12e) 

we get the elements for the Hamiltonian matrix as 

= {n + iKffjl [EphCi^d + E~ + ^{a&+ + a)&-)] \gj)\n + 1), 





= Eph{n 




-"12 












Tr{n,j) 
■"22 







(13a) 

2 V^^, (13b) 
+ 1, (13c) 
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(13d) 



Hence the Hamiltonian matrix gives 



^ /£;p/.(n + i) + ^ ^v^m 



(14) 



where I made Aj = (i?p/i — Ej) jh as the detuning of the laser from the unshifted NV center transition 
frequency. Now the eigencquation gives 

H^'^^')|vI/(-^))=i?(n,A,)|vI.(".^)), (15) 

where E{n,Aj) is the eigenvalue or eigenenergies of the system and |\I>("'-')) stands for the eigenstate or 
eigenvector in basis of IV'i" "'^) and \ip^'''^). The equation has a nontrivial solution only if 

Det(H|-^))=0, (16) 

which gives the two eigenenergies as 

E±{n, A,) =Ephin +^)± + + 1), 

^Ep,{n+\)±^^A^^+ni. (17) 

Notice that in Equ. |17[ I have defined the n-photon on-resonance optical Rabi frequency for a pulse with 
fixed duration as 

n„ = floVn + l, (18) 

so that f2„, here, has the same physics meaning as ilo in the Equs. Sll and S9 of paper which can 
be observed experimentally (thank the authors of Ref. 4J for the correction on my first version of this 
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article). Now, JIq gives the per-photon atom-photon couphng Rabi frequency. As you wiU see later, the 
experimental analysis in Ref. [4] is consistent with this coupling model. 

If we substitute the eigenenergies into the eigenequation respectively, we can solve the equation and 
normalize the results to obtain the eigenvectors finally, which depict the states of the system. We write 
the eigenstates, associated with eigenenergies E±{n, Aj), as dressed states, which read 



\n, +j) = cos($„,/2)|^("'^)) + sin(<i>„,/2)|4"^^")), 
\n,-j) = -sin(ci>„,/2)|V'i"'^'') +cos(ci>„,/2)|4"'^")), 

where j is just a label to spin states and the phase factor $„j is defined as 



— arctan 



with 



sin(i>„j/2) 



cos($„j/2) 



V2 



V2 



1 - 



A, 



A, 



1/2 



1/2 



(19a) 
(19b) 



(20) 



(21a) 



(21b) 



To get a general dynamic solution for the spin-photon system, we suppose the field is initially prepared 
in a superposition of number states 

l*/(0)> (22) 
and the electron is in excited state \ej). Thus the initial state for the whole system is 



i*(o)> = Ei*/(o))i^^) = E^»^^i")i^^-) 

j n,j 

Since Eqs. |19a| and [l9b| give 

l^^^')) =cos($„,/2)|n,+jVsin(<I>„,/2)|n,-j), 
thus the initial states for the system can be described through dressed states as 
l^'(O)) =E^".j[cos($«j/2)|n,+j) -sin($„,/2)|n,-j)]. 



(23) 



(24) 



(25) 



Because the derivation for dressed states above is in the Heisenberg picture, the states are time indepen- 
dent. If we transfer them into the Schrodinger picture, we can easily get the dynamic state vector for 
times f > as 



|*(0) = exp 



-Ht 



l*(0)), 



($„j/2)e-*^+("'^^)*/'^|n,+j) -sin($„j/2)e-'^-("^^^)*/'^|n, -j) 



(26) 



n,3 



Now let us move on to analyze the energy level shift as a characteristic of OSE. 
Equation (17 1 shows the energy splitting between the two eigenenergies are 

E,piu{n, Aj) = E+{n, A,) - E^{n, A,) = A] + Ctl. 



(27) 
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If Aj = or no detuning, the splitting energy is Espiit{n,0) = hflo\/n + 1 = firi„ (we call it the Rabi 
energy for system with n photons, i7„ is the Rabi angular frequency for this system), corresponding to the 
spin splitting between different nis states (here charaterized by flo) and associated with photon number 
state n. If Aj 7^ the energy splitting will increase, and E± moves up or down and becomes the polariton 
eigenenergy for > or Aj < to make the maximum overlap with the initial state IV'o)- 
Suppose initially the system is in ground state with total energy 

(28) 



E 1 

£^30 = Ephin + 1) - ^ = Ephin + ^) + 



After the detuning laser interacting with the NV center free-electron, the system redistributes its eigenen- 
ergies (observed energies) to E±{n, Aj). Now let's consider the case that |Aj| ^ r2„ and the photon's or 
incident light's energy is below the upper level or Aj < 0, which means the excited states will occupy the 
E-{n,Aj) level. From equation (17 1, we can get the energy shift as 



eg{n,Aj 



E-{n,A,)^ 




h\A,\ 
2 


If 




h nl 

4|Aj|' 







(A, 



(29) 



corresponding to Equ. S13 in Ref. [3]. 

Now let's consider the photon number and accumulated effects for a pumping pulse with fixed duration. 
The pulse duration r in paper ^ « I/jLs with power w l^W only contains ~ 10^ photons with wavelength 
~ 637nm and detuning in GHz range. The energy shift during the pulse lightening yields a net phase 
shift to the polariton, given by 



4A, 



where 



considering vacuum electric field 



smce 



D 



\fl\'^FDwEph cos2(5) 



\En\ = 



2nE., 



ph 



nneoAeffCT 



nEph 



(31) 



(32) 



(33) 



D/A, 



10 rad, or which is 



z rAeff 

We can estimate the accumulated phase shift for one photon is (pj 
consistent with the experimental result in Ref. j4j for Z3/27r sa lOkHz. So for the bunch of photons in 
one pulse, we can obtain an observable signal phase from the accumulated phase in the order of firad, 
as correctly analyzed in Ref. :4 by the authors. This phase change corresponds to the rotation of the 
polarized output light, so it shows the Faraday effect. 

To gain the exact expression for both the FE and the OSE, the authors used coherent state \a) to 
describe the laser field-which is valid for such laser input with considerable photon in one pulse-and 
calculated the reduced density matrices, Pspin and pught, for spin and optical components. Now suppose 
the initial state of the polariton is 

|*o)-(^/3,|5.))l«). (34) 
i 



As introduced in class, the coherent state gives 

\a) = e 



\ - a" 



(35) 
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with a = \a\e^^ and |ap — (n) bearing the mean number of photons, such that the polariton evolves to 
the state 



l*> -E/3^-^"^E47^^*^'"^^^l5.)l«)' (36a) 
— y/nl 

- E E ^-^^\9j)\n), (36b) 

= Eftl5,)|ae*^^), (36c) 

3 

where (jjj is the accumulated phase per photon by the state |(7j)|Q;). From the full density matrix of the 
resulting spin-light system, given by p = the reduced density matrix for the optical field gives 

puakt - E = E lftfl«e^*^)(«e*^^l, (37a) 

3 3 

^ \light){light\, (37b) 

with the "light" states 

-E/^J-I"^'*')' (38) 

3 

where (3j = \(3j\e''^'^i and gives the possibility of the field occupies the spin state \gj) (we will discuss 
this for detail later). 

Similarly, the reduced spin density matrix gives 

p,^,„ = {a\p\a) ^Y.f^V3 exp{-|«p[(l - e'*0 + (1 " e"'^'=)]} (39a) 

3,k 

^Y^I^kPi exp [-^{n){(|), - (l^k)] |.9,)|5fe), (39b) 

3,k 

= \spin) {spin\. (39c) 

In the last approximation, I used Taylor expansion of e"^^ and e~"^'=, since (t)j <C 1, and hence the "spin" 
states become 

\spin) = E \93) = E '^^■e^'^ • (40) 

i 3 

From the above equations, we can see that the rotated spin states are affected mainly by the detuning (A-,) 
from the light field for an given pulse and occupation possibilities on j states (associated with j3j , which 
is determined by the initially excited states using on-resonance echo technology which will be introduced 
in later sections. Since the spin probability amplitudes j3j and phase information are passed on to 
the optical field, which is described by equation ( 38 1 , through measuring the optical field's information, 
we can get the complete information on the spin-photon interaction system. This technology is named as 
"spin-light coherence for single-spin measurement and control" technology in Ref. f?]. 

Now let's consider the optical property, to build up the bridge linking between quantum optical theory 
and experimental measurement, basing on "light" states or reduced density matrix for the optical field. 
Assuming the electric field operator gives 

^(f,t) = |^ol«(r)e-^"*(a-|-at), (41) 
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where u{'r) describes the spatial mode of the optical field. The expectation value of the field gives 



{E{f,t))ught = {Ught\E{f,t)\light) 

j 



(42) 



where — + (/)^^. is the phase corresponding to the rrig = j spin state. Since phase information 
associated with cj is independent to the amplitude of the expected field value, and also a; is a very large 
number 10^^ r ad/ s), (pj dominates the observable phase information. Meanwhile, because ipi^^ only 
depends on the initial state, which can be conditioned into some fixed state, (j)^ is strongly dependent on 
= D / Aj which is basically a function of detuning Aj . And this phase information varies the magnitude 
of the measured field strength or output light density considerably (unfortunately, the authors did not 
provide this measurement data in their paper). In the far off resonant detuning experiment, only one 
linear polarization component of light is coupled to the transition channel j, hence the polarized phase is 
shifted relative to the non-interacting polarization state by an amount (pj , which can be measured by using 
an adjustable polarizing beam splitter. And the Faraday phase $i? is the difference in phase between the 
TTT-s = and TTis = — 1 (here states nig = ±1 are almost degenerate compared with rus = state) spin 
state. As the photon distribution on different polariton states is related with the electron distribution on 
different spin states (described by /3j), so that we can bring in Faraday phase ^p{N) for a system with 
total number of iV, by 

where LOg — (£'_i — Eo)/h is the frequency spacing between the resonances to iris = 0, —1 states. Assum- 
ing the high-order atom-photon interaction is negligible, the effect Faraday phase is almost a constant 
representing the per-photon effect for a system with a total photon number of TV. 

Since / oc iV, one expects to observe the FE effective phase shifting with various input light power. 
This experiment shows a good agreement with the expectation as in figure ([S]) . 

By accumulating Faraday phase in the total coherent optical field with N photons, we can get the 
relative OSE phase shift as 

<i>osE « N^F- (44) 
And the corresponding OSE frequency shift is 

^ ^OSE Pl ^ /.(.X 

= -7. = . p. ^F: (45) 

where Pl — 2nEpfi/T. So the OSE phase or frequency shift is proportional to the total laser power. 

By making the full FE line an odd Lorentzian curve, as a result of the Kramers-Kronig relation, 
and considering the dephasing mechanism is dominated by spectral diffusion and power fluctuation, the 
authors also analyzed the dephasing phenomena in experiments, which is quantized by Njj that is defined 
as the number of OSE-induced oscillations at which the amplitude envelope drops to 1/e times its value 
at T = 0. The experimental data agrees with the theoretical model quite well as shown in Fig|4] 

2.2 On-resonant dynamics theory and pulse echo technology 

Now let's specify our case in Aj — 0, that is when the light is resonant with the spinning electrons. Now 
= 7r/2, and 

E± (n, j) = E± (71, A, = 0) = Ephin + ^) ± ^^rW^ 
= Eph{n+-) ± — -, 
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Figure 3: (A)FE measurement timing sequence show the gated green laser, tunable red laser, microwave, 
and multiplexed APD timing, as well as the spin-state evolution and diode bridge signal in time. PL 
photons are binned separately in time to measure both Iple (APD measured, with Lorentz fit) and (Sz) 
(the Larmor precession rate shift due to the polarization energy shift, given by Sj = Sg{{n) /2Trh,)) for 
the rUs = — 1 and rris — prepared spin states. Gray areas indicate possible spin-polarization effects 
from the red laser, (denoted by the blue arrow at bottom-left) is the difference in measured red laser 
polarization response between the prepared = and nis — ~1 spin states. (B and C) FE data sets 
scaled identically show ^p, as well as Iple and (Sz), for both — (blue) and = —1 (red) prepared 
spin states of the ^Ey orbital- branch optical transitions at 5 and O.lmW red laser power, respectively. 
The TOs = optical transition (at OGHz laser energy) is more robust against spin polarization than the 
nig = — 1 optical transition (near SGHz laser energy), fiiad, microradians; kCps, photon kilocounts per 
second. (D) FE data set at 15mW scans across both '^Ex and "^Ey orbital-branch optical transition 
energies. The measured FE is substantially reduced for the '^Ex orbital transitions (near 16.5GHz), 
primarily as a result of the smaller optical transition energy splitting between the spin states. From 
B. Buckley, G. Fuchs, L. Bassett and D. Awschalom. Spin-Light Coherence for Single-Spin Me asurement! 
and Control in Diamond Science, 330 (2010), 1212g ]. Reprinted with permission from AAAS. 



Eq. (26) gives 



(47) 



cos(^)|V.^-))-zsin(^)|4-)) 
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Laser energy (GHz) Laser energy (GHz) 



Figure 4: (A) OSE measurement timing sequence. The red laser produces a az spin rotation proportional 
to its to 4/is pulse duration. {Sz) follows ± ^vaiaz) as a result of a ±90° microwave phase shift of the final 
Hahn echo pulse. (B) /pi of sequences that have these two microwave phases. kCts, photon kilocounts. 
(C) E5 as a function of laser power at ~ 2GHz laser energy. The error is represented by the data-point size 
and is dominated by the 15% uncertainty in the laser power calibration. The red line is a linear fit of Ss 
with 5.74MHz per microwatt slope. (D) Comparison of S5 and showing their complimentary response 
as a function of laser energy and also showing Iple and (Sz) taken with $p for both prepared spin states. 
OSE and FE data sets were taken under the same experimental conditions at 0.66-/iiy laser power. The 
Ss data point denoted by the blue arrow is the frequency fit of the data presented in Fig.3B. (E) OSE 
spin coherence measured in number (N^j) of az rotations to the 1/e decay point. Black circles represent 
the measurements with 0.66/iW red laser power (Fig. D), and green triangles represent measurements 
with 1.59/iW red laser power. The blue arrow denotes the data point fit from in Fig. B. The red line 
is a fit to the data using a model incorporating dephasing from spectral broadening and laser intensity 
fluctuations. Error bars indicate standard errors of OSE data fits. From B. Buckley, G. Fuchs, L. Bassett 
and D. Awschalom. Spin-Light Coherence foF Single-Spi n Measurement and Control in Diamond Science, 
330(2010), 1212[ 4]. Reprinted with permission from A A AS. 



According to Eq. (6a), for the n photon number state, at time t the probability of finding the system in 
1$^""''') state is 

= ^e-^'">(™+^)*/'-'C„,fcCos(%^)(ci>("'^-)|$l™^'^V™,n^. 



k.j 



g-,iJ,.(„+i)tAC'„,,cos(^) 

not. 



(48) 



= |C„,,f cos^(^). 

Since |C„j p is fixed for a given initial state of the system and can be normalized for the {n,j) states. 



the probability of measuring state is proportional to cos'^(^). Similarly, the probability of finding 



2/'f2Q^^ 
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the system in an IV'2 ) state is 

P2{n,j,t) = \Cr.jHm\^). (49) 

From the two equations above, the Rabi angular frequency flo holds for any states. The time 

evolution diagram for these two states is shown in Fig. [51 These two dynamic evolution equations (Eqs. 




Figure 5: Time evolution of a system between states iV'i"'"''') and IV'2"''''' 



48 and 49 ) form the foundations of entanglement control technology. 

For example, if at time t — to we force the system into a known state, that is C„j- is known, then 
for subsequent times the system's states will evolve with a maximum amplitude described by |C„jp. If 
we try to excite the pure = j ground state to the excited state with a resonant laser pulse of length 
T = tt/^Iq, then the system will evolve to the excited state at the end of the pulse. We call this kind of 
pulse a TT-pulse. While if the pulse length r = 27r/Oo, it has no effect on the system. So, this kind of 
27r-pulse is a transparent pulse to the system. Generally, if r is an arbitrary number, the pulse will apply 
an extra phase described by tJIq to the system. In field coupling application, the so called bang-bang 
coupling technique [43j is also developed from this phenomenon. 



3 Electron structure of NV center in diamond 

To better control the quantum state of spin and accurately realize spin-photon entanglement, we need 
to consider the fine or even hyperfine level structures of the spinning electrons of the NV center (here, 
negatively charged) in the diamond crystal environment. There have been many approaches used to obtain 
the energy band structure, for example. Local Density Function method |19j . many-body perturbation 
model ^37|, etc. In this section, I would like to calculate the matrix of the system's Hamiltonian with 
spin-spin and spin-orbit interaction, strain |17| . |39| . |53| terms, and use group theory's results to confirm 
the basis of electron states and transition selection rules with impact on the photon's polarization. 

In the absence of external strain and electric or magnetic fields, properties of the six electronic excited 
states are determined by the NV center's Csi, symmetry and spin-orbit and spin-spin interactions (shown 
in Fig|6]). Optical transitions between the ground and excited states are spin preserving, but could change 
electronic orbital angular momentum depending on the photon polarization. Two of the excited states, 
labeled \E^) and \Ey) according to their orbital symmetry, correspond to the = spin projection. 
Therefore they couple only to the |0) ground state and provide good cycling transitions, suitable for 
readout of the |0) state population through fiuorescence detection. The other four excited states are 
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entangled states of spin and orbital angular momentum. Specifically, the \A2) state has the form 

l^2)-^(|i?-)| + l) + |i?+)|-l)), (50) 

where \E±) are orbital states with angular momentum projection ±1 along the NV axis, where | ± 1) 
denotes the magnetic sublevel states with rris = ±1. Similarly, we denote mj — spin states with |0). At 
the same time, the ground states (|0), | ± 1)) are associated with the orbital state |i?o) with zero angular 
momentum projection (for simplicity, the spatial part of the wavefunction is not explicitly written). Hence, 
owing to total angular momentum conservation, the IA2) state decays with equal probability to the | — 1) 
ground state through (T-|_ polarized radiation and to | + 1) through cr_ polarized radiation (here, circular 
polarizations are represented hy a± = x ± y, while linear polarizations are represented by x and y). In 
other words, if we can trigger a photon from the electron transition from \A2) to the ground states | ± 1), 
we can entangle the frequency and spin state information into the polarized photon and easily verify the 
fidelity of entanglement by reading out the photon's polarization and frequency information and comparing 
this with our knowledge of the electron structure of the NV center. All of these properties make the IA2) 
singlet state (sometimes people use \^A2) to differentiate this singlet state from the triplet state |^A2) 
which finally forms the singlet state under interaction with the atomic environment) a good candidate for 
spin-photon entangling in diamond NV center. And, it is vital to determine the level structure of related 
states before we carry out any entanglement scheme designs. 

The inevitable presence of a small strain field, characterized by the strain splitting of {E^^y) reduces 
the NV center's symmetry and shifts the energies of the excited state (^^42)) levels according to their 
orbital wavefunctions. A group theory study of the NV center in diamond [5T], tells us that a small strain 
will not change the polarization of the transited photon, such that we can use a small strain to modify the 
energy level structure. And group theory also gives the basis of electron states and the order of eigenvalues 
as well as the transition selection rule. The basis for the electron states of interest \'^A2) gives [BD] 



\A,) = \E^)\ + l)-\E+)\-l), (51a) 

\A2) = \E^)\ + l) + \E+)\-l), (51b) 

\E^.) = l^)|0), (51c) 

\Ey) = \Y)\0), (51d) 

\E,) = \E^)\-l)-\E+)\ + l), (51e) 

\A2)^\E^)\-l) + \E+)\ + l), (51f) 



where \E±) = \ae± — e±a), \X{Y)) = \aex(y) ^ ^x{y)o), e± — — iCy and ej.{y), a span into the full 
orbitals space of the electrons. Fig. [6] shows the energy splitting for electrons in NV center by considering 
spin-orbit and spin-spin interactions. 

Including the strain effect, the total Hamiltonian reads 



El — Has + H so + H strain, (52) 

where the spin-orbit interaction term Hso, spin-spin interaction term Hgs and strain Hamiltonian are 
respectively given by jM], [BD] 

Hso = \zS,L,+\^{LS)^, (53a) 

Hss = A(|Ai)(Ai| + \A2){A2\ + m{E^\ + \E2){E2\) (53b) 

- 2A {\Ex){Ex\ + \Ey){Ey\) + A' (|A2)(^2| - \Ai){A,\) , (53c) 

Hstrain = {le^} {e^l - \ey) {ey\) + 5y {\ex){ey\ + |ej^)(e;j;|) , (53d) 



and S, L are the magnetic and orbital angular operator, the subindex _L denotes the non-axial (or x-y 
plane) component. Following the procedure as before, we can rewrite the Hamiltonian into a matrix under 
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Ai,A2, . 
El , E2 / 



A, 



Ex , Ey \ 



\ E^,Ey J2A "E. + PE, , 



Ei,E2 '*"A" 



cv-Ei - ,3-E„, 
aE2 - HE^ 



Figure 6: Diagram of energy splitting of electrons in NV center, considering spin-orbit and spin-spin inter- 
actions. From J.R. Maze, A. Gali, E. Togan, Y. Chu, A. Trifonov, E. Kaxiras, and M.D. Lukin, Properties 
of nitrogen- vacancy centers in diamond: the group theoretic approach. New Journal oj Physics, 13 (2011), 
025025 |doi:10.1088/1367-2630/13/27025025i l]. Reprinted with permission from lOP Publishing Ltd. 



the basis of [A2 Ai Ey E2 Ei]: 
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(54) 



where A = 1.42/3GHz, A' = 1.55GHz, A^ = 5.5GHz, A_l = 0.2GHz 
dence with Dr. J.R. Maze, who is one of the coauthors of reference [i 



Notice that, upon the correspon- 



Eqs. 53a together with the matrix 



in Eq. 54 should be adjusted because the transverse spin-orbit does not mix the states with different spin 
projections on the excited state. And we should use Aj_ = 7.3GHz j41) as a result of ab initio calculations 
instead of 0.2GHz which is commonly used in the early work. To perform a comparison, we maintain 
Aj^ ~ 0.2GHz and continue to calculate the eigenenergy based on the equations above. If the eigenenergy 
is A, then the eigenequation is given by 



A^ - (3 Sx^ + 3Sy^ + 64.32676667)A'^ - 85.07314251 A^ 



(3 4 



3(5„ 



6 6x^6y^ 



(82.29690576 {Sx 



123.482 {6x^ + 5y^) + 850.2580259)A2 
I- 1609.942968)A - {5x^ + 5y^) 



(55) 



3 {dy^Sx + 5y^5x ) - 59.15573333 {5x 
0.1240 5x^ + 0.3720 5x 5y^ - 787.4883241 (5^^ 



Sy'^) - 118.3114667 (5^^(5y 



740.9385799 = 0. 



I calculate the eigenenergies for states Ai, A2, Ex, Ey, E2 and Ei using Maple. By forcing 6y — 5x as 
a result of symmetry, I compared my result with Maze's result m] and another independent result [2]. 
They all agree with each other very well, except for a minor difference among the crossing points shown 
in Fig. [7] As you can see, this minor difference does not affect the rest analysis in Togan's paper, since 
only one low strain value-far less than the cross point value-is used to realize the entanglement. 

By introducing a low strain on the diamond, Togan and his teamworkers makes the energy gap between 
the 1^2) singlet state and the | ± 1) states equal to the photon energy in the laser beam with a wavelength 
of about 637.19 nm, which can be easily obtained from a commonly used tunable YAG:Nd laser. The 
level scheme used to realize spin- light entanglement is a A-type system, shown in Fig. p[a). 
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strain i' (GHz) 



Figure 7: Comparison of Electronic Structure as a function of stain. Equation [55] leads to (A), versus 
[¥T] gives (B) (From J.R. Maze, A. Gali, E. Togan, Y. Chu, A. Trifonov, E. Kaxiras, and M.D. Lukin, 
Properties of nitrogen- vacancy centers in diamond: the group theoreti c approach. New Journal of Physics^ 
13 (2011), 025025141!. Reprinted with permission from lOP Publishing Ltd.). The y-axis of Fig.(B) is 
moved down by 1.945 eV to make the and Ey originate from the origin. I have made 5y = = 0.2GHz 
in (A), which gives the similar lineshape as S± — 7.3GHz (only the slope changes slightly). The results 
of (A) and (B) agree very well, except for the crossing point under the arrow in (A). In (A), E2 starts to 
surpass Ey at the crossing point as the strain grows up. (A) much more close to the experimental result 



(yet disputable) shown in (C) (Reprinted figure with permission from A. Batalov, V. Jacques, F. Kaiser, 
IP. Siyushev, P. Neumann, L.J. Rogers, R.L. McMurtrie, N.B. Manson, F. Jelezko, and J. Wrachtrup, 
^ys. Rev. Lett. 102, pp.195506 (2009) 2 . Copyright 2009 by the American Physical Society.), where 
the solid lines are calculated with the same set of parameters in (A) , the discrete dots are measured from 
different NV centers label as NVi, NV2 and NV3. The curves in (B) do not show the energy levels' cross 
between E2 and Ey. 



By the way, the so called ^^iV and ^^C nucleus also contributes a hyperfine splitting to the electronic 
level structure in a diamond environment. To overcome this disadvantage, we can use slightly detuned 
electronic microwaves (/i waves) to cover the effect of hyperfine coupling. This fi wave usually makes a 
detuning of A = 120 MHz, and makes an OSE energy shift of fi^/A w 0.5 MHz, considering the Rabi 
frequency for this hyperfine splitting is 8 MHz. So, this effect is relatively small for the entanglement 



experiment. According to Eq.(26), by applying this jj, wave, the ground states are rotated to a rotating 
frame described by 

I ±i)t = e*'^±e-'('^±=F'^)*| ±1). (56) 
This rotation can easily transfer electrons to the |0) state from a superposition state 



|M) = 



1 

71 



1) ±e 



--i(uj-t- 



1) 



(57) 



corresponding to photon state \H) — -1=(|(t_|_) + |(t_)) and |V) = -^{\a+) — |cr_)), which are linear 
polarization states. In this way, it makes verification of non-diagram density matrix elements feasible (see 
Figj9]jb)). The supplemental materiel of Togan's paper explained this technology. 
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►Time 



Figure 8: Electron level structure and experimental procedure for entanglement generation. 

a, scheme of selected level structure used to generate entanglement: \A-i) state as the excited state and 
|±1) states as the ground states. At the beginning, after spin polarization into |0), population is transferred 
to I + 1) by a microwave 7r-pulse (i^+i)- The NV is excited to IA2) with a 637.19-nm 7r-pulse and the ZPL 
emission is collected, b, If a (t+ or cr_ photon is detected, the population in | + 1) or | — 1) is transferred 
to |0). If an \H) or \y) photon is detected, a r — 27r — r echo sequence is applied with f2+i and r2_i, 
followed by a 7r-pulse which transfers the population in \M) to |0). c, The population in |0) is measured 
using the 637.20-nm optical readout transition, d, Pulse sequence for the case where an \H) or \V) ZPL 
photon is detected (time axis not to scale). If a a± photon is detected instead, only a 7r-pulse on either 
or 0_i is used for spin readout. Inset, detection time of ZPL channel photons, showing reflection 
from diamond surface and subsequent NV emission (blue) and background counts (purple). Reprinted by 
permission from Macmillan Publishers Ltd: [Nature, 466 (2010), no. 7307, 730-734(50], copyright 2010. 



4 Decoding spin-photon entanglement experiments in NV center 
of diamond 

Now we are ready to design a spin-photon entanglement scheme for diamond NV centers. 
As discussed before, if we can form a state described by 

l*)-^(k->l+l> + k+>|-l», (58) 

where |cr±) and | ± 1) are the polarized photon and electronic spin states, then we can say we have 
formed an entangled state. In August 2010, Togan and colleagues successfully realized this entangled 
state in NV centers in diamond nanocrystal and analyzed the entanglement in a basis of four Bell states. 
The electron levels and experimental procedures used to realize this entanglement are shown in Fig. [H] 
The authors used zero-phonon-line (ZPL) photons with four basis states: cr±, = + 1^-)) 

and \y) — ^(|o'+) — lo"-)). These photon states respectively entangled with four spin states: |±) and 
|±) = ^(1 + 1) ± I - 1)). These states span the four Bell states with maximum entanglement. It is 
worth mentioning that, to read out the photon and spin states, the authors used a carefully designed 
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pulse echo sequence using the r-pulse, 7r-pulse and 27r-pulse discussed in the on-resonance model. And as 
mentioned above, they also used a temporary state \M) in a rotating frame, which can be explained using 
our off-resonant model. All of these technologies make detecting any expected states feasible. 

Mainly because of the low detection rate, they spent several weeks to do one round of experiments to 
verify this entangled state, and carefully repeated the experiments over several months to collect enough 
data to publish their results. Finally, one beautiful statics diagram is shown in Figj9] 




Figure 9: Measurement of spin-photon correlations in two bases, a, Conditional probability of 
measuring |± 1) after the detection of a (T+ or cr_ photon, b, Conditional probability of measuring |± 1) 
after the detection of an H or V photon, extracted from a fit to data shown in c and d. c, d, Measured 
conditional probability of finding the electronic spin in the state \M) after detection of a V(c) or H(d) 
photon at time t^. Blue shaded region is the 68% confidence interval for the fit (solid line) to the time 
binned data. Errors bars on data points show ±1 s.d. Combined with the data shown in a, oscillations 
with amplitude outside of the yellow regions result in fidelities greater than 0.5. The visibility of the 
measured oscillations are 0.59 ± 0.18 (c) and 0.60 ± 0.11 (d). Reprinted by permission from Macmillan 
Pubhshers Ltd: [Nature] 466 (2010), no. 7307, 730-734(60], copyright 2010. 



The probability of probing spin in | — 1) state after detecting a cr+ photon is equivalent to C{H, V) in 
Eqjlj and so on. The probability of measuring the state \M) at time t^i is 



1 ± cos(Q;(id)) 



(59) 



where a{td) — — uj-)td + (</>+ — </<-), with fixed 4>± as the initial condition for every experiment. We 
can see the results give a very high fidelity for verifying the entanglement states. 



5 Conclusion and outlook 

In the approaches presented in this paper, we find the Jaynes-Cummings model can be successfully applied 
to analyze entanglement dynamics and technologies, including phase rotating, energy splitting (using strain 
or optical pump) and pulse echo technologies, in diamond NV centers, which is a promising medium for 
practical Quantum Informational and Quantum Computational applications. Our analysis of the electron 
structure also works well in this case. These approaches and technologies can be potentially used in similar 
materials and systems. 

Through out the discussion, we can see that a well selected and controlled level structure is the key 
to realizing spin-photon entanglement. To make the spin information map onto photons, it is necessary 
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to choose two spin opposite states and a common ground or excited state, which is the so-cahed A-type 
structure. In fact, some other level structures are also discussed in recent publications for realizing such 
an entanglement. For instance, a M^-type structure is clearly discussed here [28]. To fully understand this 
level design technology, a tomography analysis may be useful. 

At the same time, to improve the entanglement performance, coupling the NV center mode to a high-Q 
cavity is also helpful (note from June, 2011: just reported by Patton and O'Brien that this enhancement 
was realized in a diamond microring [46]), especially to improve the entanglement distance and hence 
realize multi-qubit entanglement in solid-state materials. And since realizing entanglement networks is 
the final goal, it is necessary to study the coupling effects based on cavity-QED theories and discover 
practical quantum devices based on entangled units. To realize this goal, the symmetry of multiparticle 
entanglement [53] and many-body dynamics theories [37] may also be helpful. 
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